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THE DIRICHLET PROBLEM FOR THE CONVEX ENVELOPE 



ADAM M. OBERMAN AND LUIS SILVESTRE 



Abstract. The Convex Envelope of a given function was recently charactor- 
^ ^ ' ized as the solution of a fully nonlinear Partial Differential Equation (PDE). 

Cn ' 111 this article we study a modified problem: the Dirichlet problem for the 

underlying PDE. The main result is an optimal regularity result. Differentia- 
^J 1 bility (C^'" regularity) of the boundary data implies the corresponding result 

for the solution in the interior, despite the fact that the solution need not be 
|/~\ ' continuous up to the boundary. Secondary results are the characterization of 

the convex envelope as: (i) the value function of a stochastic control problem, 

and (ii) the optimal undercstimator for a class of nonlinear elliptic PDEs. 

< 
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1. Introduction 



In this article we study a problem which is at the interface of Convex Analysis 
and nonlinear elliptic Partial Differential Equations. The object of this study is 
the Convex Envelope. In a previous work by one of the authors |16| . a Partial 
^ ' Differential Equation (PDE) in the form of an obstacle problem for the convex 

CO i envelope was obtained. In this article we further explore the connection between 

C^ ' the convex envelope and nonlinear elliptic PDEs by studying the Dirichlet Problem 

for the convex envelope. 

The convex envelope, u, of given boundary data, g, on a bounded domain, 
f^ I fl C M", with boundary, d^, is the solution of the following fully nonlinear and 

degenerate elliptic Partial Differential Equation, 



(PDE) AiM(x) = 0, for a; in rj, 

(D) u{x) — g{x), for x e dil. 



Here Ai[u] denotes the smallest eigenvalue of the Hessian, D^u, of the function u. 
j^ ' The solution describes a convex function which is nowhere strictly convex, in a 

sense which we explain later. We refer to the solution of (|PDE[) . (|D|) as the convex 
envelope of the boundary data .g(x), since these solutions agree with the usual 
definition (|CE|) (below). 

While the convex envelope has been well-studied, it is useful to have a different 
characterization. In fact, this characterization has led to PDE-bascd numerical 
methods for computing the convex envelope |17j . 

Continuity at the boundary. Note that while we study the Dirichlet boundary condi- 
tion ini for (|PDEp , the solution need not achieve the boundary values (for example 
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if the boundary data is nonconvex). In addition, in dimensions three or higher, the 
convex envelope itself need not be continuous up to the boundary |14j . 

Interior regularity. The first regularity result establishes C^'" regularity of the so- 
lution in the interior, provided the boundary data is C^'". 

This result is unusual for the following reason: loosely speaking, regularity re- 
sults for elliptic PDEs (see [5] or [15] for more references) often fall into one of 
two categories. The first is regularity from the boundary. For example, solutions 
of uniformly elliptic PDEs arc C^'" regularity up to the boundary, provided the 
boundary data is C^'" (Theorem 9.31 of [U]). 

The second is interior regularity, assuming nothing on the boundary nevertheless 
forces some regularity in the interior. Our first regularity result is of the first type. 
However it has the unusual feature that despite the smooth boundary condition, 
the regularity holds only in the interior. 

Geometry of the contact set. The second result we obtain concerns the geometry 
of the contact set of the solution with a supporting hyperplane. As a consequence 
of this result, it is established that at each point in the interior, there is at least 
one direction in which the solution is linear. This result is used to show that 
viscosity solutions satisfy (jFDEp in a classical sense even though they may not be 
differentiable. 

Stochastic control interpretation. Wc give an interpretation of the convex envelope 
as the value function of a stochastic control problem. 

Application to degenerate elliptic PDEs. Pursuing the connection with elliptic Par- 
tial Differential Equations further we show that the convex envelope of Dirichlct 
data is a natural subsolution of a class of degenerate elliptic PDEs. 

Viscosity solutions of (|PDE|) arc also viscosity solutions of the elliptic Mongc- 
Ampcrc equation with zero right hand side. Thus the results for (jPDEp also apply 
for the elliptic Mongc- Ampere equation with zero right hand side. 

The equation (|PDE[) is highly degenerate: it is degenerate in every direction 
except the direction of the eigenvector of the first eigenvalue of the Hessian. Another 
highly degenerate equation is the Infinity Laplace equation [l] , which is degenerate 
in every direction except the direction of the gradient. 

1.1. The Convex Envelope. The convex envelope of a given function g{x) is 
mathematical object which has been the subject of study for many years. The 
convex envelope of the function g(x) is defined as the supremum of all convex 
functions which are majorized by g, 

u{x) = sup{w(a;) | v convex, v{y) < g{y) for all y e M"}. 

It was recently observed [TB] that the convex envelope is the solution of a partial 
differential equation. 

The equation for the convex envelope, u, of the function g : M" -^ M, is 

(1) max{u(a;) — (/(x), — Ai[u](a;)} = 0. 

Here Ai[u](a;) is the smallest eigenvalue of the Hessian D'^u{x). See Figure [T] 
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Figure 1 . Illustration of the equation 



The equation ([T]) is a combination of a fully nonlinear second order PDE, Ai [u] {x) , 
and an obstacle term. In dimension n = 1, this equation coincides with the classi- 
cal obstacle problem. But the differential operator is fully nonlinear in dimensions 
n = 2 and higher and it is degenerate in all but one direction. 



1.2. Related results. The possibility of an equation for the convex envelope was 
suggested by [T^]. A computational method for computing the convex envelope 
which used a related equation was perfomed in [20]. Methods for enforcing con- 
vexity constraints in variational problems have been studied as well [6|. For more 
references on computational work, see [17j . 

The regularity of the convex envelope u{x) of a given function g{x) has been 
studied by [131 [T^ [5]. In this case, it has been shown that when the envelope 
function g{x) is C^'", the convex envelope is as well. The analysis for the envelope 
problem is somewhat easier, since supporting hyperplanes will touch the function 
g at some point where it is differentiable with matching derivatives. 

In dimensions three or higher, the convex envelope of a function defined on a 
closed convex set need not be continuous up to the boundary [M] . 
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2. A Partial Differential Equation for the Convex Envelope 

The convex envelope can be defined on a domain fl in M" even when g{x) is 
defined only on part of the domain. This is achieved simply by setting g to be 
infinity outside the domain of definition. 

Although much of our arguments still hold for closed sets C, there is the possi- 
bility that the envelope is infinite in parts of the domain. Since here our interest 
is in the relationship of the convex envelope with elliptic PDEs, we restrict our 
attention to case where C is 9r2, the boundary of fl, and assume g{x) is defined 
and bounded on C. 

Then we define the convex envelope of the boundary data g{x) to be 

(CE) u{x) = sup{v{x) I V convex, v{y) < g{y) for all y S dVi] 

In this last definition, the set of convex functions v can be replaced with the set of 
affine functions, since a convex function is the supremum of its supporting planes. 

The function defined in (|CE[) will not usually be continuous up to the boundary, 
since g{x) may be concave at some points. However, we cannot simply take a con- 
vexification of g and expect the same result, since g is defined on lower dimensional 
set. 

We now focus our attention on the Dirichlet problem (jPDEp .([D |) . The prob- 
lem (JPDE|) . ([D|) can be obtained from ([T]) by setting g{x) to be infinite on the 
interior of the domain. 

Analogous results for concave functions can be made by replacing Ai with A„, 
the largest eigenvalue. 

2.1. Convexity. For basic definitions of convexity, see [4] or the appendix of [9]. 
For more on convex analysis see the textbooks [3| or [19] . The function u : R" -^ M 
is convex if for all x,y £ i?" and < i < 1 

u{tx + (1 - t)y) < tu{x) + (1 - t)u{y). 

If u : M" — >■ M is convex, then for each xq G K" there exists a supporting hyperplane 
to u at xq. In other words, there exists an affine function P{x) = u{xo)+p- {x — xq) 
such that 

ju{x) > P{x) for all x e R" 

\u{xo) = P{xo) 

If u is differentiable at x, then Vu(a;o) = VP(a;o) and P is unique; if not, there 
may be more than one supporting hyperplane. 

For twice-differentiable functions, convexity can be characterized by the local 
condition that the Hessian of the function be nonnegative definite, 

D^u{x) > for all x £ M" if and only if u is convex. 

Note that the first condition is equivalent to Ai[u](a;) > 0. The characterization 
is valid even for continuous functions, provided the equation is interpreted in the 
viscosity sense, as will be explained in the next section. 

2.2. Viscosity Solutions. The theory of viscosity solutions is a powerful tool for 
proving existence, uniqueness and stability results for fully nonlinear elliptic equa- 
tions. The standard reference is the User's Guide [8]. A readable introduction 
is the Primer, by Crandall [7]. The theory applies to scalar equations of the form 
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F[u] = F{D^u{x), Du(x), u(x), x), which are nonlinear and ehiptic, i.e. nondecreas- 
ing in the first argument. Solutions are stable in the sense that if the equations 
F"^ converge to F, the corresponding solutions u'^ converge to the solution u of 
F, uniformly on compact subsets. Uniqueness is a consequence of the comparison 
principle: if F[u\ > F[v] in D,, with u > v on dCl, then u > v in CI. Viscosity 
solutions can be constructed using Perron's method, 

u{x) = sup{0(a;) | is a subsolution of F[(p] = 0}. 

This last construction will coincide in our case with the definition of the convex 
envelope (|CE|) . The definition of viscosity solutions for (jPDEp follows. (In the 
interest of clarity, we omit a standard technical argument which requires working 
with upper and lower semi-continuous envelopes, see Section 4 in [8].) 

Definition 2.1. The continuous function u is a viscosity subsolution of (jFDEp if 

u{x) < g{x), for x e dfl 

and if for every twice- differentiable function (t>{x), 

(2) — Xi[(j)\{x) < 0, when x is a local max of u — (j). 
The continuous function u is a viscosity supersolution of (jPDEp if 

u{x) > g(x), for x G dQ. 
and if for every twice- differentiable function (j){x), 

(3) — Ai[(/)](x) > 0, when x is a local min of u~ (j). 

The function u is a viscosity solution of (jFDEj) if it is both a subsolution and a 
supersolution. 

The following theorem was proved in jl7) . 

Theorem 2.2. The continuous function u : R" — > R is convex if and only if it is 
a viscosity subsolution of (jPDEp . 

An immediate consequence of the previous theorem and Perron's method is the 
characterization of the convex envelope of the boundary data as a viscosity solution 
of (jPDE|) . This result was given in more detail in [16] for the convex envelope of a 
function (rather than the Dirichlet data). 

Compare the Perron formula 

u{x) = sup{w(a;) | v{x) is a continuous subsolution of ()PDEp }. 

and the definition (jCEp . Use the definition of viscosity solutions and Theorem 12.21 
to see that the supremum is over the same set of functions. (Here semi-continuity 
of the function u(x) defined above follows from the definition. Continuity can be 
established by an additional technical argument, see Theorem 4.1 in [5].) 

2.3. Application to Monge-Ampere equation with zero right hand side. 

The Monge-Ampere equation with zero right hand side is 

det{D^u) = 

In order for the equation to be elliptic, there is additional constraint that u be 
convex. While the next result is probably not new, the characterization of the 
convex envelope (jPDEp gives a clear and concise proof and statement. 



6 ADAM M. OBERMAN AND LUIS SILVESTRE 

Lemma 1. Viscosity solutions of the elliptic Monge- Ampere equation with zero 
right hand side are also viscosity solutions of (jPDEp 

Proof. Using Ai [u] > to enforce the convexity constraint, and writing 

Xi[u] < ■ ■ ■ < Xnu[u] 
for the eigenvalues of the Hessian of u and noting that 

det{D^u) = Xi[u] ■ ...Xn[u] 

we see that 

det(Z?^u) = 0, and u convex 
is equivalent to the simpler condition 

Ai[it]=0. D 

3. Geometry of the Contact Set 

The geometrical characterization of the convex envelope ([T|) suggests that solu- 
tions of (|PDEp are convex, but nowhere strictly convex. By studying the geometry 
of solutions we can give a rigorous meaning to this statement. 

It is easy to build examples to show that solutions of (|PDE[) need not be differ- 
entiable; for example u(x, y) = \x\ is a solution. Neither must they be continuous 
up to the boundary. For example, when fJ is a square in the plane, centered at the 
origin and g{x, y) = x^ — y^, the solution of (|PDEp is u{x, y) = x^ — 1. 

For classical solutions, by which we mean twice continuously differentiable (C^), 
the principal curvatures of the graph of the solution must be 0, A2, with A2 > 0. 
This means that at every point in the domain, solutions are flat in the direction of 
the first eigenvector of the Hessian. 

However for nonclassical solutions, it is possible that this direction is changing. 
The direction of flatness has implications for the regularity of solutions and for the 
stochastic control interpretation. 

Nonclassical example solutions have the property that at any point x there is 
a direction d = d{x) on which the solution restricted to the line i{s) — x + ds \s 
linear. 

Is it true that for all points x in the domain, there is a line segment 
containing x on which the solution is linear^ 
To answer this question, we investigate the geometry of the contact set. 

Definition 3.1 (Contact Set). Let u he a solution of (|PDEp and let P^ he a 
supporting plane (jSHp to u at x. Define the contact set to he the set of points where 
the supporting plane touches the graph of the function. 

(4) C^^{ye^\u(y) = P'^{y)}. 

The contact set may also depend on the choice of the plane P^, if there is more 
than one supporting hyperplane at x. 

3.1. Statement and Proof of the Theorem. 

Theorem 3.2. Let ft he a compact domain in M" and g : dfl -^ M. Let u : ft ^ M. 
he its convex envelope. For any point x inside fl, let P^ he its supporting plane and 
C^ its contact set. Then 
(1) C^ is convex. 
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(2) C^ intersects the boundary offl. 

(3) C^ contains a line segment connecting x to a point on the boundary of fl. 

(4) C" is the convex hull of C" n dVt. 

Proof. (1) For any y, z € C^ , let w ^ (y + z)/2. Since u is convex, 

^^ u{y) + u{z) _ P-{y)+P-{z) 

On the other hand, since w S C^, and P^ is a supporting hyperplane, 

/ X P''(y)+P'''(z) 
u{w) > P{w) = ^-^^ —. 

So C^ is convex. 

(2) Suppose C^ does not intersect the boundary of il. Since C^ is convex, it is 
closed, so there is a small neighborhood D of C^ which also does not intersect the 
boundary of il. Since P^ is a supporting hyperplane, u > P^ in I? \ C^. Let 

e = min u{y) — P^{y). 
y&dD 

Then e > 0. So define u{y) = max{u{y) , P^ (y) + e/2). Then u{y) > u{y) at some 
point y in D \ C^ , but u{y) = u{y) outside D. This contradicts the definition of 
u, (|CE| . 

(3) Suppose C^ intersects the boundary of il at y. Then since C^ is convex, 
there must be a line segment from x to y. 

(4) Let D be the convex hull of C^ H dil. Since C^ is closed, so is D. Moreover, 
since C^ is convex, D C C^. We are left to prove the opposite inclusion. 

Let y £ il\D. We will show that y ^ C^ . Since C^ is closed and convex, then by 
the Plane Separation Theorem [3] there is an affine function L such that L{y) > 
and L < in D. 

Since L < in D, then g — P^ is strictly positive in the compact set 9r2n{L > 0}. 
Therefore, there is a 5 > such that g - P"^ > S in dfin {L > 0}. Thus, for e > 
small enough, P^ + eL < g on dil. 

Since u is above any affine function that is below g on 5il, then u{y) > P^{y) + 
sL{y) > P'-^{y). So y ^ C^ , which finishes the proof. D 

Let X £ C^ , then by the last part of the theorem, 

(5) a; G conv[yi,...,yfe], ioi y^ e dft 

for i = 1, . . . , fc with 2 < k < n (where conv[yi, . . . , yk] means the convex hull of 
yi). This means that u is affine on a (A; — 1) dimensional set, which contains x in 
its relative interior. 



3.2. Consequences of the Geometry. Theorem 13.21 means that, even if the 
solution is not differentiable, it always satisfies (|PDEp in an (almost) classical sense. 
Recall the definition of the second directional derivative. 

Definition 3.3 (Second directional derivative). Given a unit vector \v\ = 1, the 
second derivative of u in the direction v is given by 

d^u h f u{x + vk) — u{x)\ k f u{x) ~ uix ~ vKV 

lim 111 



dv'^ h,k^o+ h + k\ k J h + k 
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Corollary 3.4. Let u{x) be the convex envelope of the boundary data, defined on 
a domain fl. Then 

d'^u, 



Xi\u\ix) ~ mill -—^ix) = 
|'u|=i dv^ 



for all points x inside Vi. 



Proof. For all ft,, k, using the definition of a second directional derivative as a limit, 
followed by the minimum characterization of the smallest eigenvalue of the Hessian, 
gives: 

h f u(x + vk) ~ u(x)\ k f uix) — u(x — vhY 

\i[u\{x) = mm hm ' i i i 



\d\=ih,k^Q+ h + k \ k J h + k 

Since C^ is the convex hull of C^ n dfl, u{x) is linear on a fc — 1 dimensional afhne 
set containing x, which is given by ([5|). So choosing any direction d which lies in 
this set, we get zero in the last equation. D 

4. Stochastic Control Interpretation 

In |16] the convex envelope was reinterpreted as the value function of a stochastic 
control problem. For the Dirichlet problem (jPDE[) . an even simpler interpretation 
is available. Our derivation is formal but can be made rigorous; we refer to |10| for 
a rigorous derivation of related equations. For readers not familiar with stochas- 
tic control problems, an introduction to optimal control and viscosity solutions of 
Hamilton- Jacobi equations can be found in [9]. 

4.1. The Control Interpretation. Consider the controlled diffusion 

(d^it)^V2e{t)dw{t), 

\ x(0)=a-o. 

where w is a one-dimensional Brownian motion, and the control, O(-), is a mapping 
into unit vectors in R". The process stops when it reaches the boundary of the 
domain, at which point we incur a cost g{xr)^ where r is the time when it reaches 
the boundary. The objective is to minimize the expected cost 

J{x,0{-))=E-[gixr)] 
over the choice of control 9{-). The value function is 

u(x) = min J(x,6'(-)), 

0(.) 

which describes the minimal expected cost at a given initial point x, assuming that 
an optimal control strategy was pursued. 

In the general setting of stochastic control problems, the value function satisfies 
a fully nonlinear Partial Differential Equations, which is obtained by applying the 
Dynamic Programming Principle (DPP). The DPP provides a link between nearby 
values of the value function by assuming a constant control pursued over an infini- 
tesimal time interval. A readable introduction to this principle in the deterministic 
case can be found in [9]. 

Now apply the DPP to derive the equation for the value function. One strategy 
is to fix 9{-) = 6* to be constant, and let the diffusion proceed for time t, thereafter 



THE DIRICHLET PROBLEM FOR THE CONVEX ENVELOPE 9 

following the optimal strategy. This strategy costs E^°'^[u{x{s))] = E^°[u{x{s)) \ 
9{-) = 9]. Minimizing over 9 gives 

u{xa) == min^^«^^[w(x(t))] + o{t). 

Q 

Using the definition of infinitesimal generator corresponding to the diffusion ((S]), 
with 9 fixed (see e.g. [H]), gives 

(7) lim^l>WM^^.,T^V.),. 

Using ([7]) in the preceding equation and taking the limit t — > yields 

- inf 9'^ D^u(x)9 ^0, 

"1=1 

along with u = g on the boundary. Finally, using the Raleigh-Ritz characterization 
of the eigenvalues 

Xi{M) = min x^ Mx, 

\x\ — l 

for symmetric matrices, M, recovers (jPDEI) . 

5. C^'°' Regularity 

In this section we prove a regularity result. We begin with an example to show 
that the interior regularity result is optimal in the sense that it cannot be continued 
up to the boundary. 

We will use the notation 

Br^{xe R" I |xp < r}, Sr^{x(=, R" I Ixp = r}. 

for the ball of radius r, and the sphere of radius r, respectively. 

5.1. Optimal interior regularity example. Let Q, = Bi d M^, the unit ball in 
two dimensions. Consider the function 

u{x,y) = -{l+xf-' 

for e > small. The function is convex and continuous on Bi. Writing 

u{x,y) ^ f{r,9) = -{l + rcos9f~^ 

in polar coordinates, we see that u restricted to the unit circle is C^'^~^'^ , since near 
9 = zbvr the function behaves like 9^~^'^ . Compute the derivative, 

^/(I,0) = (l-e)sin0(l+cos0)-^ 

which has a singularity about 9 = ±7r which is on the order of 9^^^'^ . 

Notice that the function -^u{x,y) has a singularity at (—1,0). So the function 
u(x, y) is not C^'" up to the boundary of Bi for any a > 0. 

Thus we have provided the desired counterexample to the regularity of the so- 
lution of (|PDEp . (|D|) near the boundary. 
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5.2. Proof of the Regularity result. Let n = Bi C M" and let g : Si C M" -^ E 
be a C^'" function. Consider (|PDEl) .(P) 

{Ai[w](a;) = 0, for x in Bi 
u = g, for X on dBi = Si 

This equation is understood in the viscosity sense. Reinterpreting the definition 
from the Perron construction gives the characterization 

{u{x) = max {L{x) : L{x) = A-x + b,L<g on dBi} , for x in Bi 
u ^ g for X on dBi = 5*1 

We recaU that if the boundary data g is the restriction of a nonconvex function, 
then the solutions will not be continuous up to the boundary. 
We next prove the interior regularity result. 

Theorem 5.1. Consider (jPDEp . ((D|l . for Q = Bi C M". If the boundary data g is 
C^'" on dfl ~ Si for some a G (0, 1], then the convex envelope u is C^'" in -B1/2. 
Moreover, the following estimate holds, 

ll"llci.°{Si/2) - '^ll5llci-°(9Bi) 

where the constant C depends only on the dimension n and on a. 
Proof. 

Preliminaries, Resettling. We can assume 

llffllci.-OBi) = 1 

because if we can prove the estimate in this case, the general case will follow by 
rescaling. 

To prove the estimate, we will choose two points, xi,X2 G Bii2i £^iid establish 
the estimate 

\Vu{xi)-yu{x2)\ <C\xi~X2r 

for some constant C depending only on the dimension n and a. First replace the 
gradients with the gradients of supporting hyperplanes. So let 

(8) L^ = A^■{x-x^) + b,, i = l,2 

be the supporting hyperplanes to Xi, respectively. The argument will show that the 
supporting hyperplanes are unique. Define 

M{x) = ma.x{Li{x),L2{x)). 

See [Figure 2[ a). 

By subtracting an appropriate affine function from u and g we can assume with- 
out loss of generality that 

Ai+A2 = 

and 

(9) min M = 0. 

x£Bi 

so the {n — 1) dimensional hyperplane P = {x \ Li{x) = L2{x)} can be written as 

P = {x\ M{x) = 0} 
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M = max{L,, y 




Figure 2. (a) The function M = max(Li,L2) (b) P separates 
x\^xi and meets S\ at a bounded angle. 



and note for later that we can write 
(10) 



Mix) = ^— ^dist(2;,P) 



As a result of these simplifications, our goal is to show that 

(11) |Ai-^2|<C|xi-X2r 

for some constant C depending only on the dimension n and a. 

The proof will proceed based on an estimate on min g from below in terms of 
1^1 — ^2 1, which will lead to an estimate on \x\ — X2\ from below. 

Comparing Euclidean and geodesic distances. We will be comparing distances from 
points on the boundary S'l = dBi to the intersection of P with 5*1, both in the 
ambient space M" and along the boundary, using the geodesic distance. With this 
in mind, write dist(j/,.T) for the standard Euclidean distance, and write for y € Si, 

dists,(y,PnS'i) 

for the geodesic distance between y and P n 5*1 on 5'i . 

Since P separates xi and X2, P passes through S1/2 and intersects Si nontan- 
gentially at some minimum angle Oq > (which depends only on the dimension n) . 
Refer to [Figure 2[ b). Thus for any point y G Si, the distance between y and P is 
greater than, but comparable with the geodesic distance between y and P D dBi 
on dBi. We record this assertion as follows. There is a constant C2 > 1 such that 
for any y G Si 



(12) 

where 

(13) 



dist(2/,P) < dists,(y,Pn5i) < C2dist(y,P) 



1<C. 
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Setup estimates using distances on the boundary. Choose xo G Si so that the mm- 
imuni of g is achieved at xo , 

(14) min g{x) = g{xo) > 

and moreover, this minimuni is nonncgative, by ^. 

Now select a great circle on 5*1 which is perpendicular to P and which passes 
through xq. Let j{t) be a parameterization by arclength of the great circle, which 
passes through P and xq 

7(0) €P, l{ti)=xo. 

Using (1121) . noting the parameterization is by arclength, we obtain 

(15) dist(7(i),P) <i < C2dist(7(i),P), for < i < 7r/2 

where we have chosen 7(0) in the same hemisphere as xo- 

By convexity of u and the supporting hyperplane conditions, 

g(x) > u{x) > M{x), for x in Bi. 

For future reference, define the function of one variable G{t) = ,9(7(i)) and record 
this last result as 

(16) G{t) > M{j{t)), for < t < 7r/2 

Apply estimates at two points. We will proceed to apply (jlSp along with (jl6|) at ti 
and at a second, larger value of t to obtain the desired result. 
Combining ([T5)) with the expression ([TU]) for M{x), we obtain 

(17) dist(xo,P)< ^^^^"^ 



\Ai-. 

Applying (fT5|) at ti we obtain 
(18) ti<C2dist{xo,P) 

Combining (fT8|) and (|17p we obtain 

2g(a:^o) 



(19) O^^^^^^A,-^, 
Next let 

(20) ^. = 2C,,'^^"") 



1^1 - A2 



then 

(21) C,^^<t.-t,<2C,,'^("") 



1^1-^1 - 1^1 
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Then compute 




Gih) > Mhit2)) 


from ([TOl) 


-l^^;^^ldist(7(t,),P) 


using (liO|) 


^ \Ai-A,\ 1 ^ 
- 2 Cj' 


by (|15| 


IA1-A2I 1 ^^ 2g{xo) 
2 C,'\A,^A,\ 


inserting (|20|) 


= 25(2:0) 


after simplifying 


which we record as 




(22) Gih) > 2g(xo) 





Conclusions from the estimates. Using ([22]) and the definition G{ti) = 5(2:0) we 
obtain 

'''G'(t)dt = G(t2)-G(ti)>5(xo) 

Next we use the fact that J^^ f{t) dt > /3 imphes there exists t* in [ti,<2] such 
that f{t*) > /3{t2 — ti)~^. Apply this fact to the previous integral to conclude the 
following. There exists t* in [<i,t2] such that 



(23) 



G'{t*)>gixo)it2-hr'> ^ '^^ ^'' 



2G2 



where we have again used (PT|) . 



Relate to assum,ptions on 5. Since 5(0:0) is the minimum of 5. 

G'(ii)=0. 
Moreover, since Ii5||pi,<:,(a5 ) = 1 then 

|G'(si) - G'(s2)| < \si - 521", for ah si,S2 
In particular, 

G'{t*) < ir-tiT 

so substituting ([25)1 into the last equation 

' '^^-^<r-t,r<|t2-t,r<f2G2,'^^"") ^" 



2G. 



1^1-^2! 



using (121]). This last equation simplifies to 5(^0) > (4G2)"^+^/"|^i - ^2^+^/^ 
Thus, there is a constant c such that 



min5 = 5(.2:o)>c|Ai-A2r+i/^ 
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Conclude. Therefore the constant function L{x) = c\Ai — A2\^^^^°' is below g on 
5"!, and then, by comparison 

u{x) > c\Ai - ^sT+^Z" for X in Bi 

Returning to the original points xi and X2, and recalling the supporting hyper- 
planes dl]), we have 

u{x,) = M{x^) = -L^— -^dist(a:,, P), for i = 1,2 

then 

dist(x„P)>c|Ai-^2r/", 
which clearly implies (|lip which was the desired result. D 

6. Estimators for elliptic PDEs. 

6.1. Sharp underestimators for elliptic PDEs. Convex functions are natu- 
ral candidate subsolutions for nonlinear elliptic partial differential equations. The 
equation (|PDE|) allows us to prove that the convex envelope is the sharpest subso- 
lution for a class of equations, in the sense we describe below. 

Similar results about best overestimators can be obtained for the concave enve- 
lope of the boundary data, which is the solution of — A„[u] = 0. 

What we describe below is closely related to the maximal and minimal Pucci 
operators, which are the best possible sub- and super-solutions for a class of uni- 
formly elliptic equations with given ellipticity constants, as in [S]. We begin by 
showing that the convex and concave envelopes provide the best possible sub- and 
super-solutions for a class of possibly degenerate elliptic equations. We then review 
the uniformly elliptic case, following [5]. Finally we show that the Pucci maximal 
and minimal operators converge to the convex and concave envelope as the ratio of 
constants goes to zero or infinity, respectively. 

To that end, consider the Dirichlet problem for a nonlinear elliptic partial dif- 
ferential equation, P[M](a::) = F{D^u{x), Du{x),u{x),x). For motivation, suppose 
that we have very little detailed information about the operator F. Assume for 
clarity that F satisfies 

(24) F{Q,p,r,x) = Q, for all p € M",r G ]R,x e f7. 

Also assume that we have (precisely known or estimated) Dirichlet boundary con- 
ditions, 

u = g, on dQ. 

Then we show below that we can estimate the solution by the convex and concave 
envelope of the boundary data. 

The case ([M|) can be generalized to F{0,p,r,x) = /, with corresponding results 
hold for the operators \i[u] = f. 

6.2. Estimators without ellipticity bounds - envelope operators. Without 
assuming that the operator is uniformly elliptic, we obtain the equation for the best 
underestimator. As expected, wc get the convex envelope operator. 

Define the best underestimator for elliptic equations 

(25) u = inf{M(x) I F[u] = in fi, u = .g on dn, F satisfies (p4)) }. 

F 
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The boundary conditions hold in the viscosity sense, which is why the result may 
not achieve the boundary conditions at all points. 

Theorem 6.1. Consider the class of nonlinear elliptic equations which are homoge- 
neous in the sense of ()24p and satisfy the boundary data g. The best underestimator 
for this class is the convex envelope of the boundary data. 

Proof. We make note of the fact that the convex envelope of the boundary data, 
which we write as u'^, is the solution of (jPDEp with boundary data g. 

Now simply apply the definition ([25]) with F = —Xi to obtain u<u'^. 

To obtain the other inequality, first verify that the convex envelope is viscosity 
subsolution of the equation. We do this by checking Definition 12.11 Suppose x is 
a local max of u'^ — cj), for some C^ function cj). By Definition 12. 1[ Ai[(?!)](x) > 0, so 
D'^(I>{x) > 0. Now compute 

F[4>]{x)=F{D^(l}{x),D4>{x),(l){x),x) >F{0,D(l){x),4>{x),x) =0, 



which follows since F is elliptic, and by ([24]). Now since u'^ is the convex envelope 
of the boundary data, u'^ < g on dVt. Together these facts imply that u'^ is a 
subsolution of F . The comparison theorem applied to F yields u'^ < u. 

Together, these results imply u = u'^. D 

6.3. Estimators with elHpticity bounds - Pucci operators. We next state 
the case where we have some additional information on the operator F . Suppose 
in addition that 

(26) F[m] is uniformly elliptic, with constants < 7 < nV. 

This section follows [5], which should be referenced for the definition of uniform 
ellipticity in the fully nonlinear case, and for additional details. (For the fully 
nonlinear case where F is differentiable, uniform ellipticity is equivalent to the fact 
that the linearization of the nonlinear operator F at any particular values should 
have eigenvalues bounded above and below by the ellipticity constants.) 
Define the best underestimator for uniformly elliptic equations, 



i7,r 



= mi{u{x) I F[u] = in O, M = 5 on dVt, F satisfies ^ and ([26l)}. 



Then there is an explicit equation for u_, p . It is the solution of the Pucci minimal 
equation, [5l sec. 2.2], which can be written as 



^u. 



(27) M-^M = -(7^A, +r^A,], {AOr=i eigenvalues of D 

In the two dimensional case, the result is the operator 

■^7,rN = -(7-^1 +rA2), Ai < A2 eigenvalues of D^u, f7 C M^. 

This follows from the simple observation that in order for the A^r^p[u] = 0, the 
eigenvalues must have different signs. 

There is a related equation for the best overestimator, which is the Pucci maximal 
equation, see sec. 2.2]. 



16 ADAM M. OBERMAN AND LUIS SILVESTRE 

6.4. Convergence of Pucci Operators to the convex envelope. We are now 

equipped to prove the following. 

Theorem 6.2. The best underestimator for uniformly elliptic equations converges 
uniformly on compact subsets to the best underestimator for elliptic equations, which 
is the convex envelope of the boundary data, as T/j — >■ oo . 

Proof. The result follows from stability of viscosity solutions, once we show that 
■^7 rl"] ^" ^^i["] ^^ r/7 -^ 00. Divide through by T in (f27|) . and take the limit 
r/7 -^ 00, to give X^A <o '^« ^ ^' This last identity is equivalent to miuj Ai = 0, 
which recovers Ai = 0. D 
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